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A classical analysis suggests that an external magnetic field can cause trajectories of charge carriers 
on a superconducting domain wall or cosmic string to bend, thus expelling charge carriers with 
energy above the mass threshold into the bulk. We study this process by solving the Dirac equation 
for a fermion of mass m/ and charge e, in the background of a domain wall and a magnetic field of 
strength B. We find that the modes of the charge carriers get shifted into the bulk, in agreement 
with classical expectations. However the dispersion relation for the zero modes changes dramatically 
— instead of the usual linear dispersion relation, ujk = k, the new dispersion relation is well fit by 
u! ~ TO/ tanh(fc/fc*) where fc* = to/ for a thin wall in the weak field limit, and fc* = eBw for a 
thick wall of width w. This result shows that the energy of the charge carriers on the domain wall 
remains below the threshold for expulsion even in the presence of an external magnetic field. If 
charge carriers are expelled due to an additional perturbation, they are most likely to be ejected at 
the threshold energy ~ to/. 



I. INTRODUCTION 

Fermion zero modes on domain walls and cosmic 
strings have been known to exist since the 60's 0,12. If 
the fermions are electrically charged, it is generally pre- 
sumed that the string or wall is superconducting in the 
sense that it can sustain long-lived persistent currents 
(3 . Such superconducting strings could be present in the 
universe, and would have observable electromagnetic and 
other signatures (for a review, see Q). Light supercon- 
ducting strings could also be trapped in the center of 
our galaxy and could explain Q the observed positron 
annihilation line 

The observable signatures of cosmic superconducting 
strings depend on the magnitude of the currents they 
carry. As the strings move under their tension, they 
cut across ambient astrophysical magnetic fields and this 
generates currents on the strings. Interactions of charge 
carriers on the string are expected to limit the maximum 
current that a string can carry 01 ■ Hence it is important 
to understand the interaction of a superconducting string 
with an external magnetic field. 

In this paper we consider a straight, static, current- 
carrying, string, placed in a background magnetic field. 
The magnetic field is assumed to be perpendicular to 
the string. The charge carriers on the string feel the 
background magnetic field via the Lorentz force which 
is also perpendicular to the string. Thus the ambient 
magnetic field tends to bend the charge carriers off the 
string as depicted in Fig. ^ 

The dispersion relation for zero modes on a string (with 
no external magnetic field) is that for a masslcss particle 

LOk = +k (1) 

where, to be specific, we have taken the -I- sign. (If the 
zero mode propagates in the opposite direction along the 
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FIG. 1: The charge carriers run along the string. In the pres- 
ence of a perpendicular external magnetic field, the Lorentz 
force would tend to bend the charges off the string. 

string, we would have uj^ = — fc.) The dispersion relation 
for the remaining ("non-zero") modes is of the form 

LUk = ±Vk2 -hTO2 (2) 

where m is a mass associated with the mode. For scat- 
tering states, m is the mass of the fermion in vacuum, 
m f . For a particle in a zero mode state to escape, it has 
to transition into a scattering state. Since a particle in 
a zero mode only has kinetic energy, to escape into the 
bulk, its kinetic energy has to transform into mass energy. 
For k > ruf, this is not forbidden by energy conservation 
because it is possible to find the scattered momentum, 
k', such that = + mj. 

Let us consider a particle in a zero mode with k/m 
oo. The particle certainly has enough energy to escape 
the string. The question is whether an external mag- 
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netic field can trigger this escape. At first sigiit it would 
appear that the fermion zero mode is very fragile since 
it is a bound state as regards its motion perpendicular 
to the string, yet it is degenerate with continuum states 
that are unbound. Such bound states arc known to dis- 
solve resonantly into the the continuum under the effect 
of even a very weak perturbation, for example, in the 
auto-ionization of excited atomic Helium. However, as 
shown by Fano 0, in order to destroy a bound state 
a perturbation must couple it to continuum states with 
which it is degenerate. Due to the translational symme- 
try of the problem it is easy to see that a uniform mag- 
netic field does not couple the zero mode to continuum 
states that are degenerate with it^. One way to under- 
stand this result is to note that all the energy of the zero 
mode is kinetic energy, while only a fraction of the parti- 
cle energy off the string is kinetic. To be ejected from the 
string, some of the kinetic energy of the zero mode has 
to get converted into mass energy, and the magnitude of 
the momentum of the particle off the string therefore has 
to be less than that of the zero mode. Thus we are lead 
to the conclusion that at least perturbatively a uniform 
magnetic field does not destroy the zero mode by Fano's 
mechanism. It is therefore necessary to analyse the ef- 
fect of the magnetic field non-perturbatively in order to 
definitively settle the question. In this paper we carry 
out such a non-perturbative analysis. 

For the reader who may not wish to go through all the 
details, we briefly summarize our result for the dispersion 
relation for fermions on a string or domain wall in the 
presence of a magnetic field. The dispersion relation for 
small wavevcctors, fc, resembles that of a zero mode on 
a string without any magnetic field. Intuitively, k is the 
momentum of the particle and at zero momentum, there 
is no Lorentz force and so the magnetic field plays no role. 
Hence the dispersion relation is linear for small k. On the 
other hand, at large fc, the magnetic field dominates and 
the dispersion relation is as if there is no domain wall. In 
this regime, the energy is independent of fc, just as in the 
Landau problem. The transition of the dispersion from 
linear in k to independent of k occurs when the effect 
of the domain wall is comparable to that of the external 
magnetic field, or approximately at fc* where k„ = mf 
for a thin wall in the weak field limit, and fc* = eBw for 
a thick wall of width w. Thus the fermions in the lowest 
energy level have a dispersion relation that is of the form 

ujk = TO/tanh^^^ (3) 



^ The situation is different if the external magnetic field is 
non-uniform. However, in an astrophysical context, the non- 
uniformity occurs on astrophysical scales and the ejection rate is 
correspondingly small. Other effects that can eject zero modes, 
such as string curvature, also occur on these scales. Here we will 
only be interested in whether there are effects that can cause 
faster ejection than that given by astrophysical scales. 



This dispersion is exact in the thick wall limit and fits 
our numerical results reasonably well in other cases. 

The modified form of the dispersion relation for su- 
perconducting strings may be important in astrophysical 
applications as briefly described in the concluding sec- 
tion. 



II. SIMPLIFYING THE PROBLEM 

The proper framework for discussing the effect of ex- 
ternal magnetic fields on zero modes on strings is in terms 
of the field theory of scalar and gauge fields that gives a 
string solution, the fermions which have zero modes on 
the string, and the external magnetic field. Such a model 
is quite complicated and it is desirable to find a simpler 
model with all the relevant physics. For this reason we 
consider a field theory in two spatial dimensions that has 
domain wall solutions, a fermion with zero modes on the 
domain wall, and a magnetic field. The Lagrangian for 
the model containing a scalar field, (j)^ a fermion, ^j, and 
a vector field, A^, is 

+ i^-/''idf, + ieAf,)ij- gcjj^iP (4) 

The 7 matrices are defined, in terms of the Pauli spin 
matrices, by 

7* = , 7^^ = io-2 , 7^ = -icn (5) 

The scalar sector of the field theory is the A(/)^ model 
and has the classical domain wall solution 

(/)(a;) = jytanh ( — ] (6) 

where rj = m/VA and w = \/2/m is the width of the 
wall. 

In the absence of an external magnetic field, the 
fermionic zero mode is 




This solution has zero momentum and zero energy. It 
can be boosted along the direction of the wall to get a 
propagating zero mode 

i>k{t,x,y) = cxp [-i{u!t - ky)] ^|:o{x) (8) 

The dispersion relation is as in Eq. |^ 

col = +k (9) 

In Sec. lIVI we will discuss the zero modes when we also 
have an external magnetic field. To gain intuition for 
the final results, we now discuss the Landau problem of 
fermions in a magnetic field but without a domain wall. 
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III. RELATIVISTIC LANDAU PROBLEM 

We now consider the fermion modes in a magnetic field 
taken to be uniform in the xj/— plane. The potential is 



= (0,0, -Ba;) , n = t,x,y 



(10) 



With this choice of gauge, the fermion wavefunctions 
can be factorized as in Eq. (jSJ and we can solve for the 
functions tp^ . 



There is a special eigenmode with l^Jqi^ 



(11) 



where 



Xo = 



eB 



(12) 



The remaining relativistic Landau eigenfunctions are 

j^/2eB(n+l) 



A. 



\^ mf±y/mj+2eB{n+l) 

y^Q-i{'^'i!.,k,±*-^V) 

with energy eigenvalues 



LPn+l {x + Xo) 



=±^m2 +2eB(n + l), n-0,1,2,. 



(13) 



(14) 



Note that the energy eigenvalues in the Landau problem 
are independent of the wavevector k. 

Here An^± are normalization constants and Lpn{x) is 
the solution to the shifted harmonic oscillator equation 
which can be expressed in terms of Hermitc polynomials: 



Vn{x) 



{eB) 



1/4 



7ri/42"/2n! 



Hn ( xy eB ) exp 



-eBx' 



(15) 



IV. DOMAIN WALL IN A MAGNETIC FIELD 

Now we consider the full setting where there is a do- 
main wall along the y— axis, and a uniform magnetic field 
of magnitude B in the a;y— plane (see Fig. EJ. Then the 
scalar field is as given in Eq. © and the gauge field in 
Eq. 

The Dirac equation now reads 



{il^d^, + eBx-^^ - g(t>{x)) V' = 



(16) 



This is the equation for a Dirac particle moving in a 
uniform magnetic field with a position dependent mass 
m{x) = g(j){x). Far from the wall, the mass approaches 
the asymptotic value ±m/ as a; ^ ±oo, where mj = grj 
To find the eigenstates of Eq. H16(l , we write 



(17) 



Domain wall 




FIG. 2: In our simplified setup, a domain wall lies along the 
?/— axis and a uniform magnetic field lies in the xj/— plane. 
Fermion zero modes run along the wall i. e. in the j/— direction. 



The equation for the x fields is 



■m{x) — to 
—idx + ik + ieBx 



—idx — ik — ieBx 
~m{x) — LO 



Xi 

X2 



(18) 

The dimensionful scales that characterize the system are: 
w, nif, and eB, with w and m/ both occurring within 
m{x). These dimensionful parameters can be written 
as three length scales in the problem: w, l^, = "mj^ 

and Ib = 1/VeB. In addition, a mode is characterized 
by its wavelength 27r/fc. Only the ratio of the various 
length scales is important. For convenience, we work in 
units where w = 1 except in the thin wall approximation 

0. 

This 



(Sec. IIV D|) where we take w 

We would like to find all solutions to Eq. H18|) . 
is possible to do analytically in two extreme limits (1) 
thick domain wall (w 3> lip, Ib), and (2) thin domain wall 
{w <^ l^, Ib)- In general, the modes need to be found 
numerically. Before discussing the solutions, it is neces- 
sary to clarify the concept of a zero mode in the presence 
of a magnetic field and to prove that such a mode does 
exist in general and not only in the two limiting cases for 
which we can find exact solutions. 



A. Zero mode at fc = 

At zero magnetic field, the key distinguishing features 
of the zero mode state are that it is localised in a direction 
perpendicular to the domain wall (along the x-axis) and 
has zero energy at /c = 0. With a magnetic field turned 
on, in the present gauge, all wavefunctions are localised 
along the x-direction, but we can still ask whether there 
is a zero energy state and, if so, identify it as the mag- 
netic field analog of the zero-mode state. To this end, we 
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must study Eq. H18() and ask whether it has well-behaved 
solutions with lu = when we set fc = 0. 

In the zero field case studied in section^ we were able 
to establish the existence of the zero mode by explicit, 
exact solution of the Dirac equation. This is a rare lux- 
ury; under more typical circumstances it is necessary to 
resort to index theorems Q or other indirect arguments 
to establish the existence of zero modes. The strategy we 
follow here is to eliminate X2 and rewrite Eq. H18() as a 
second-order equation 



92 



to' d 



ox"^ TO ox 



w? + eB + e^B^x'^ 



eBx XI = 0, (19) 



where to = m{x) and to' = dm/dx. Once this equation 
has been solved for xi we may determine X2 via 



i d eB 

X2 = -^Xi — xxi 

m ox TO 



(20) 



For X —^ 00 we may set m' = in Eq. 119() and 
thereby find that in this asymptotic region xi obeys 
the Schrodingcr equation of a non-rclativistic oscillator. 
Therefore it has one solution that grows as .x 00 and 
one that is well-behaved and decays. The key question 
now is whether a well-behaved solution can be found at 
a; = that will connect with the well-behaved solution at 
-|-oo. Since to(x) = for x = 0, it is clear that the origin 
is a (regular) singular point for the differential equation 
(I19|l and it appears to be a possibility that of the two 
independent solutions near the origin only one is regular. 
However, by use of a series expansion, we can show that 
there are in fact two independent solutions that are both 
regular at a; = 0. One has the series expansion 



Xi = l + lU-2^eB- 
8 V Ml 

the other has the expansion 



e^B^ 



Xi 



M3 4 , 

—X + 
2/ii 



(21) 



(22) 



Here we have taken m{x) = IJ.1X + ji^x^ + . . . near the ori- 
gin. A generic combination of these two solutions would 
develop into a superposition of the growing and decaying 
solutions as X — > 00; but by fine-tuning the coefficients 
we can arrange for the solution to decay as x ^ 00. Fur- 
thermore, since the solutions in Eqs. (|21|l and H22() are 
both even, and the differential Eq. 1)19(1 is also even, it 
follows that the solution will be even and well-behaved 
as .X —> —00 as well. This concludes our demonstration 
that there is a zero energy solution to Eq. H18() for fc = 0. 



B. Zero mode for non-zero k 

We turn now to the magnetic zero mode for non-zero fc. 
It is shown below that for a fixed fc Eq. H18|) has a discrete 



spectrum of eigenvalues to. Thus we may hope to identify 
the zero mode at non-zero fc as the state that smoothly 
evolves into the zero energy state as fc 0. That such 
a state exists is not obvious: there is a competing possi- 
bility that the spectrum is not continuous as fc — > and 
that the zero energy state at fc = has no counterpart 
for non-zero fc. In this section, however, we will give a 
variational argument that makes it certain beyond rea- 
sonable doubt that the spectrum is continuous at fc = 
and that the zero mode therefore continues to exist for 
non-zero fc. Further, we are able to verify this continuity 
and to study the evolution of the zero mode with fc ex- 
plicitly in the thick and thin wall limits discussed in the 
subsections below. 

To study the magnetic zero mode for non-zero fc it is 
helpful to rewrite Eq. (|18|l as a second-order equation by 
analogy to Eq. H19|) . We obtain 



dx^ 



d 



{lj + to) dx 



Xi + W - ni^ 



{uj + to) 



(fc + eBx) 



eB - {k + eBxf 
Xi = (23) 



where we have used 

i d 



ieBx 



X2 



Xi 



{uj + to) dx {lj + to) 



■Xi+ikxi (24) 



to eliminate X2 from the coupled Eq. (|18|l . As x ±00 
we may set m! — Q and find that Eq. (|23() turns into the 
equation of a non-relativistic oscillator. Therefore it has 
a discrete spectrum. 

For Lo outside the mass gap (w > to^ or a; < —rrif) Eq. 
(|23|l is not singular on the real axis. Thus for any uj we 
can find a solution that is well behaved at the origin and 
as X — > 00. By tuning w to a quantised energy level we 
can also ensure that it is well behaved as x — )■ — 00. Thus 
it is easy to see that a discrete spectrum will arise outside 
the mass gap. However for the present purpose these 
states are of no interest since they will not continuously 
evolve into the zero energy state as fc — > 0. 

If we look for states that lie below the mass gap 
(—TO J < LO < ruf) we must come to grips with the fact 
that Eq. (|23|) is singular at xq where xp is fixed by the 
condition u + to(xo) = 0. By making a series expansion 
about this point it is easy to show that Eq. (|23|l has one 
well behaved solution at this point and one that is sin- 
gular. Thus we must fine-tune w in order to obtain a 
solution that behaves well as x 00'^. Since Eq. (|23|l is 



^ More precisely, we find that both solutions for xi are regular at 
XQ, but when we compute X2 using Eq. 1241 we find that one 
of them leads to a singular X2- By tuning w we can ensure that 
both solutions at xg are regular and hence that the solution at 
X ^ 00 is well behaved. Alternatively we can accept that only 
one solution is regular at xq. For generic lu this solution will 
correspond to a superposition of growing and decaying solutions 



5 



not symmetric about xq, it is not obvious that the solu- 
tion will be well behaved as x — *■ —oo. This raises the 
real concern that there may be no eigenstates within the 
mass gap for non-zero k and therefore no continuation of 
the fc = state of zero energy. 

To put this concern to rest and to establish that there 
arc in fact states within the mass gap we turn to a vari- 
ational argument. To this end, it is convenient to write 
Eq. lfTH|) in the form 

Here H is an hermitean matrix differential operator and 
we regard /c as a fixed parameter. Now has a spec- 
trum that is bounded from below and it is therefore an 
operator to which the variational principle may be ap- 
plied. We will seek a normalised state jfi) that satisfies 
{n\H'^\n) < mj. This implies that H'^ has an eigenvalue 
below mj: for every value of k, and in turn that H has 
an eigenvalue within the mass gap^. 
For the state we take 

A straightforward computation reveals 

-XrifllN'^in) = 3-41n2- — fln2- i I 
mj nif V 2 J 

+ J_(l + e + le'BA. (27) 
mj \ 2 7 

We would like to ensure that the right hand side of 
Eq. H27|) is less than unity. For suitable parameters this 
is clearly the case. Thus at least over a certain range of 
TO/ and eB there are sub-gap states for small non-zero fc. 
By continuity we may therefore reasonably conclude that 
quite generally there is a state that is the continuation 
to non-zero fc of the zero energy state at fc = 0. We will 
confirm this expectation by approximate solution to the 
problem in the thick wall limit and exact solution in the 
thin wall limit below. 



C. Thick wall 

In the limit that the wall thickness is large it is possi- 
ble to derive asymptotic expressions for the zero mode as 
well as the lowest few excited states that lie outside the 



as x — ♦ oo; but by tuning ui we may expect to be able to arrange 
for the solution regular at xq to evolve into one that is well 
behaved at a; — » oo. 
^ Since [H, H-^] = 0, eigenstates of can be chosen to be eigen- 
states of H and the eigenvalues of are simply the squared 
eigenvalues of H. 



mass gap. To understand the approximation involved it 
is instructive to set to' — > in Eq. (|^ . On doing so 
we see that xi obeys the Schrodinger equation for a har- 
monic oscillator centered at xq = —k/eB. The width 
of the oscillator wave function is 1/ \J eB. Since we arc 
working in units where the wall width it; = 1. it would 
appear that the approximation of treating m{x) as con- 
stant and equal to m{x^ should be accurate when the 
extent of the wavefunctions is small compared to the dis- 
tance over which the mass varies, i.e., in the high field 
regime XjieB) ^ 1. In Appendix fXl we will show that 
the condition for validity is slightly more restrictive. 

In summary the thick wall approximation consists of 
replacing m{x) in Eq. H18|) with the constant to(xo) where 
Xq = —k/{eB). The solution to Eq. H18|l in this approx- 
imation may then be read off from the solution to the 
relativistic Landau problem given in section 11111 Pro- 
ceeding in this way, we find that the zero mode state 
that lies within the mass gap has an energy 



uo%{k) = —m{xo) ~ mf tanh ( — ] (28) 



and a wavefunction given by Eqs. Hll() and 115|1 . 

There are several noteworthy features of this zero- 
mode solution. First observe that indeed uig as 
fc 0, showing that this state is continuous with the 
zero energy state at fc = 0. Hence this state is identifi- 
able as the zero mode state for non-zero fc. Next we see 
that in contrast to the situation at zero magnetic field, 
the dispersion relation of the zero mode is linear only for 
small fc. For larger fc, saturates and always remains 
within the mass gap even as fc ^ ±oo. At zero field the 
zero mode wave function is always localised about the 
domain wall; indeed it is independent of fc in this case. 
By contrast with a magnetic field the zero mode wave- 
function is centered about xq = —k/{eB) and it moves 
far away from the domain wall as fc ^ ±oo. The shape 
of the wavefunction is always a Gaussian and does not 
depend on fc. 

Evidently in the thick wall limit the low lying excited 
states of Eq. ifT^ are also given by Eqs. 1^3} , and l(T^ 
if we substitute to — > m{xo) = — to/ tanh(fc/ei?). For 
sufficiently excited states this approximation will break 
down because the width of the wavefunctions grows to a 
point that it is no longer appropriate to treat m{x) as a 
constant equal to its value at the center of the wavefunc- 
tion. 

To conclude we give the conditions for validity of the 
thick wall approximation to the zero mode. It is shown 
in Appendix 1X1 that the approximate solution is good if 
eB :S> I and either |xo| = \k/eB\ ^ 1 or to/ <^ eB. 
Here, as elsewhere in this subsection, we have worked in 
units where the wall width w = 1. 
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D. Thin wall 

111 this approximation the waU profile corresponds to a 
sharp step 



m{x) ~ mfQ(x) — mfQ{-~x) 



(29) 



With this profile Eq. (|23|l leads to the following second 
order equation for xi 



_ m) -eB- [k + eBxf] xi = 



x^Q 
(30) 



This is the equation of a shifted harmonic oscillator, 
h the 
= kit 



With the change of variables z = \/2eB {x + xq), where 
xq = kJeB, Eq. H3()() takes the form of Weber's equa- 
tion 



1 



z 

T 



XI =0 



Here 



2eB 



(31) 



(32) 



is not generally an integer. 

A well-behaved, i.e. tending to as x — > ±oo, solution 
to Eq. H31|) can be given in terms of the Weber function 



D 



y {^/2^{x + XQ) 



CDy {-V2^{x + xo) 
where v is given in Eq. (|32|) and 



a; > 
a; < 



(33) 



C 



Dy {xoV2^ 



Dy (-XoV2eB 



(34) 



is found by demanding continuity of the wavefunction at 
the origin. Note that the derivative dxXi will be dis- 
continuous at X = due to the non smooth wall profile 
there. The quantization of in Eq. H32|) arises in the 
simple harmonic oscillator problem in quantum mechan- 
ics from imposing that dxXi be continuous, but in the 
problem at hand, the energy will only be determined af- 
ter requiring continuity of the lower component, X2, of 
the eigenmodes 



X2{x) = 



uj-\-m f 



[dxXi - (fc + eBx)xi] , x>0 
[dxXi - (fc + eBx)xi] , X <0 



(35) 



X2{x) = 



-iC 



^/^iBvDy^i - 2{k + eBx)Dy 



^/2lBvDy^i - 2{k + eBx)Dy 



, X < 
(36) 




FIG. 3: Dispersion relation for the zero modes and the lowest 
few excited states that lie outside the mass gap for eB — ruf = 
1. The integers n± correspond to the integer n in Eq. (I14II 
and the ± sign that occur in the relativistic Landau problem. 
Note the symmetry (k) = — a;„_ {—k), which is exact and 
not limited to the thin wall approximation. 



where C is also given by Eq. (|34|) . 

Requiring that X2 be continuous leads to the quantiza- 
tion relation that determines the possible energies of the 
eigenmodes: 



AkrufDy [-ky/2/eB] [ky/2/eB] = -V2^ 



X \^{lu + m)Dy^i i^-k^/2/eBj [k^/2/eBj 

+ {uj- m)Dy_i (kyj2/eB^ [-ky/2/eB^ (37) 

where v is given in Eq. (|32|l . 

In order to numerically solve the condition in Eq. 137|l , 
it is convenient to express the Weber functions in terms of 
hypergeometric functions. This is done in Appendix IbI 
In Fig. 121 we show the shape of the dispersion relation 
for the zero modes for eB ^ nif = 1 and a few of the 
scattering states. 

The numerical results from Eq. 1)37(1 can be approxi- 
mately fit by Eq. Indeed, using the results in Ap- 
pendix we can study the dispersion relation around 
the origin, fc = 0, to find fc*: 



k^ 



dk 



^l — v 



k=0 



This can also be expressed in terms of Weber functions 

x>0 



(i + „2)r(i + i^) 
1 



r 1 



(38) 



where v = nif /v2eB. The absolute value of the slope 
du / dk at the origin is always less than 1, approaching 
this value asymptotically for large values of i.e. for 
small magnetic fields. 
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V. CONCLUSIONS 

We have studied fermionic zero modes on a domain 
wall with an external magnetic field in the configuration 
shown in Fig. [21 This is more tractable than the case 
of zero modes on a string in an external magnetic field 
(Fig. ^1 and is expected to have similar physics since the 
essential ingredients of the two systems are identical. 

We have found that the zero mode continues to survive 
even in an external magnetic field but is now centered off 
the wall. This is perhaps not so surprising because even 
at the classical level the magnetic force tends to deflect 
the fermion away from the wall. On the other hand, we 
are accustomed to thinking of the zero mode as being 
localized on the wall and, in the external magnetic field, 
this is certainly not the case. 

In the absence of an external magnetic field, the dis- 
persion relation for the zero modes is linear and the en- 
ergy increases indefinitely with increasing momentum. In 
the presence of the external magnetic field, this picture 
changes dramatically and the dispersion relation has the 
shape shown in Fig. The form of the dispersion rela- 
tion for the zero mode can be approximated by 



1 , 

ojk = TO/tanh — 



(39) 



where, in the case of a thick wall, Eq. H28|) shows that 
this approximation is exact with fc, = eBw, w being the 
wall thickness. For the thin wall, the fit is more accurate 
for small values of w, as shown in Fig.01 where it is also 
shown that fc* can be approximated by its exact value at 
the origin given by Eq. Note that |wfc| < m for all 

k. It is also worth reiterating that the dispersion relation 
is linear only for k < k^,. For larger k, the dispersion sat- 
urates and hence charge carriers actually propagate with 
a lower group velocity, dtok/dk, which vanishes exponen- 
tially fast with increasing k. In the astrophysical context 
magnetic fields are weak in the sense nif ^ V2eB and 
so the crossover scale fc* « m,/. 

For a fermion in a zero mode to escape the wall, it 
necessarily requires some additional energy to transition 
from the zero mode to an energy larger than mf. The 
additional energy can be exponentially small for k > k^ 
and hence these fermion zero modes are extremely fragile 
- even a small perturbation can knock out the particles in 
these modes. For weak perturbations, the particles will 
be ejected into the lowest non-zero (n 7^ 0) modes which 
also have energy ^ mf. The perturbation itself can be 
inhomogeneities in the magnetic field, string motion and 
curvature, or the interactions of particles outside and on 
the string. 

The modified dispersion has consequences for the cos- 
mology of superconducting strings, including the forma- 
tion of stable vortons In connection with the sce- 
nario of positron production by superconducting strings 
in the Milky Way [3 , our analysis indicates that positrons 
will be produced at an energy comparable to their rest 
mass, or about 511 keV. Then the annihilation of these 




FIG. 4: The slope, dw/dfc|^^Q = m//fc*, as a function of 
V = m//\/2e_B. The symbols show mf/k, as found by a 
global fit to Eq. 13911 for ruf — 0.1, 0.5 and 1 (circles, triangles 
and squares). The curve is a plot of Eq. and there is 
good agreement with the global fit for small v and for all m/ . 
The deviation at large v may be due to numerical issues and 
departures from the fitting formula in Eq. (I39II . 



positrons with ambient electrons will also produce 511 
keV 7— rays. 
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APPENDIX A: VALIDITY OF THE THICK WALL 
APPROXIMATION 

We now evaluate the leading corrections to the thick 
wall approximation used in section IVC. To this end it 
is helpful to write the Hamiltonian in Eq. (|18|l in the 
form H = Hq + V. Here Hq is the Hamiltonian with the 
replacement m{x) 'm{xo) where xq = —k/{eB). The 
perturbation 



V = 



m{x) — m{xo) 

^{xo) — m(x) 



may be further decomposed as V 
Vi = m'{xo){x - xo) 



V2 



1 



m"{xo){x - xqY 



Vi + 


1 

-1 



(Al) 
where 

(A2) 



In the thick wall limit we expect that m(x) is slowly 
varying so that V2 is nominally smaller than Vi . However 



8 



we will find that both perturbations have a comparable 
effect on the energy of the zero mode state. 

The zeroth-order solution to H is the thick wall ap- 
proximation presented in section IVC. The perturbation 
Vi does not perturb the energy of the zero mode to first 
order since (zero |Vi| zero) = 0. A straightforward calcu- 
lation shows that the second order correction due to Vi 



zero 



^ tanh xosech'^xo . 



(A3) 



To assess the accuracy of the zeroth-order solution, we 
note that 

ruf — ujg 2mfCxp{—2xo) , a;o — *■ cxd (A4) 
On the other hand the correction 



,(2) 



3 

-nij exp 



oo 



(A5) 



Thus the perturbation is far smaller than the distance 
to the mass gap and may be considered negligible for 
Xo 3> 1. We would also like to ensure that the pertur- 

(2) 

bation Wzoro does not shift the zero mode energy outside 
the mass gap for xq ~ 1 or Xq <C 1. Therefore our crite- 

(2) 

rion becomes |u;zoro| <^mf, which leads to the condition 
ruf <^ eB. If Vi was the only perturbation we would 
conclude that the zeroth-order approximation is accurate 
when Xo 3> 1 or when niQ eB. 

However, we must also take into account the perturba- 
tion V2. Although nominally smaller than Vi , V2 perturbs 
the zero mode energy to first order and therefore has a 
more significant effect. The first order correction is 



^zoro = -^CXp(-2.To). 



(A6) 



Applying the same reasoning as for the second order cor- 
rection we conclude that the perturbation V2 is insignifi- 
cant when either (i) xo 3> 1 and eB ^ 1 or (ii) eB ^ ruf 
and eB ^ 1. These conditions subsume the conditions 
for Vi to be negligible and are therefore the conditions 
that must be satisfied for the validity of the thick wall 
approximation to the zero mode. 



APPENDIX B: FROM WEBER TO 
HYPERGEOMETRIC FUNCTION 

We can find the dispersion relation of the zero modes 
for a given magnetic field by finding the roots of Eq. I|37() 



that satisfy w < m/ for different values of k. For this pur- 
pose it is useful to consider the expansion of the Weber 
function in the complex plane in terms of the Kummer 
confluent hypergeometric function iFi{a; p; z): 




1^ 1 2' 2' 2 



2/ K 2) 

This representation is valid along the entire real axis. 
Since it is simpler than any form we have found in the 
literature, we briefly sketch its derivation. We start by 
observing that the Weber Eq. may be transformed 
to the confluent hypergeometric equation by changing to 
the independent variable s = 2:2/2 and the dependent 
variable ^ = e'^^^xi- By means of this transformation, 
which is valid for z > 0, we obtain a confluent hyperge- 
ometric equation [H p. 671] with parameters a = —1^/2 
and b = 1/2. Thus we obtain the solution to Eq. 



1/ 1 Z2 

2' 2' T 



(B2) 



Here U2 is a confluent hypergeometric function of the 
third kind [Hp. 612], which is well-behaved as z — > 
00, and the constant pre-factors have been inserted to 
be consistent with the usual convention m p. 1641]. 
We may now use a well known joining formula |lll p. 
673] relating U2 to the standard confluent hypergeometric 
function iFi to obtain the final result HH1(I . 

So far, our derivation applies only to z > 0, but we 
can see from the differential Eq. (|31|l that Di,{z) is ana- 
lytic at z = 0, and, since the function iFi is known to be 
analytic over the entire complex plane, we conclude that 
Eq. (|Bip holds along the entire real axis. We note that 
the standard treatises 0, both arrive at the repre- 
sentation Eq. Hijl() . valid for z > 0, but offer a much more 
complicated continuation of it along the negative z-axis 
(see, for example, the asymptotic analysis in sec- 
tions 16.51 and 16.52, or the identical discussion by 
on p. 1641). 
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